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ABSTRACT 


It is noted that certain common linear wave operators have the 
property that linear variation of the initial data gives rise to one- 
dimensional evolution in a plane defined by time and some direction 
in space. The analysis is given for operators arising in acoustics, 
electromagnetics, elastodynamics, and an abstract system. 
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1 Introduction 


In this paper we point out an interesting fact concerning three common instances of linear wave 
propagation (acoustic, electromagnetic, and elastic). In each case, we consider the initial-value 
problem in some region of space sufficiently small that the gradients in the data can be taken 
as constant. Then the evolution of that data in time corresponds precisely to that produced by 
one-dimensional waves travelling in some distinguished direction. In the three instances, this 
direction is defined 1 * by the pressure gradient, by (V x B) x (V x D) (not the Poynting vector), 
and by any vector perpendicular to V x v. We also consider an abstract system, to which all 
2x2 systems in a certain class are equivalent. 

2 Statement of Result 

Consider a set of linear partial differential equations in n unknowns and d space dimensions, 

i=d 

d t q + = 0. (1) 

i=i 

A simple wave solution of ( 1 ) is a solution of the special form 

q(x, t) = ref(l-x- At) 

where l is a unit vector in the direction of wave propagation, A is the speed of propagation, and 
/ is an arbitrary scalar function. Such solutions exist if and only if r/ is a right eigenvector of 

t=d 
1 = 1 

Linearly varying initial data for (1) can be considered as a vector space V in R nd , conveniently 
displayed as d column vectors 

D = [d Xl q,d X2 q,-] 

Data giving rise to a simple wave is of the special form 

Dl = * * ’]• 

Denote by Si the subspace of V containing such data, with i fixed and A / 0, i.e., the set of 
data giving rise to one-dimensional wave motion in the direction i . 

The evolution operator A that produces dtq from an element of V is an n x nd matrix whose 
nullspace Af{A) comprises the set of linearly varying steady solutions to (1). We claim that for 
some interesting evolution operators, such as those mentioned in the Introduction, an arbitrary 
element of V can be represented as a sum of elements lying in Af and elements lying in some <S*, 

£ = + ( 2 ) 

j k 

1 The notation used in this paper is that bold type represents a vector such as a set of unknowns; vectors with 

geometric meaning in R 3 have arrows surmounting regular type 
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Note that we do not claim that AT + Si = D, but that (when the trick works) the direction £ 
can be chosen so that A f + «S/ includes any given element of D. To evaluate the coefficients in 
(2) we operate on both sides with A, giving 

&t<l = ^2 <*k^kTt,k (3) 

k 

which is merely a decomposition of the time derivative onto the eigenvectors for wave motion 
along t. This differs from the one-dimensional decomposition (Riemann problem) in that the 
direction of i is unknown and enters into the equations nonlinearly. In the following sections we 
see how this works out for various special cases. 


3 Examples 

3.1 Examples with Unique Solutions 

3.1.1 Acoustics 


For the equations governing acoustic waves in a uniform medium we have 


At = 


0 l\ (2 1 3 

0 0 0 

i 2 0 0 0 ’ 

4000 


( 4 ) 


where the unknowns are (p,v) T , the pressure and velocities, and units have been chosen to 
make the sound speed unity. The eigenvalues are Ai i2 , 3 ,4 = -1,0,0, 1 and the non-stationary 
eigenvectors are 



1 


1 

ri = 

-t 

•-i 

II 

e 


The expansion (3) is 


which can only be solved if 
when the wavestrengths are 


<*1 


Pt 

v t 




Ml = — gradp, 


2 (t'Vt-Pt), 


«4 = - $ + Pt)- 


( 5 ) 


3.1.2 Electromagnetics 


Maxwells equations in free space, for the magnetic field 

B , and the electric field D , in units for 

which the speed of light is unity, give rise to 







■ 0 

0 

0 

0 

—^3 

^2 



0 

0 

0 

^3 

0 



A,= 

0 

0 

0 

^3 

0 

—12 

-t 2 
0 

0 

0 

0 

(6) 



0 

1 1 

0 

0 

0 



ll 

-lx 

0 

0 

0 

0 
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The eigenvalues Ai, 2 , 3 , 4 , 5,6 are -1, -1,0,0, 1, 1 and the non-stationary eigenvectors are 



s 


t 


t 


s 

r l = 

t 

, r 2 = 

s 

» r 5 - 

—s 

, r 6 = 

-t 


where s,t are any pair of unit vectors such that £, s, t form a right-handed orthogonal system. 
The expansion (3) is 


implying 


B t 

Dt 

= CL 1 

—s 

-t 

+ a 2 

-t 

—s 

+ a 5 

t 

—s 

+ 06 

s 

-t 


(?) 


B t = (a 6 - 01)5 + (<15 - a 2 )t 

D t = -(c*6 + )t- (as + <* 2)5 


For a solution to exist, both Bt and Dt must lie in the plane spanned by s and t , which is the 
plane normal to t. Then 

1 || B t x D t = curl# x curl/5 (8) 


and we easily find 


Ql = 

— Bt ■ s — Dt ■ t, 


CL2 = 

-B t • t- D t ■ s, 


^5 = 

Bft — Df $, 

(9) 

<*6 = 

BfS — Dft. 



3.2 Examples with Nonunique Solutions 
3.2.1 ‘Cauchy-Riemann’ Equations 

Consider the system obtained by adding time-dependent terms to the Cauchy-Riemann equa- 
tions; 


U t + U X + Vy = 0 

v t - V X + Uy = 0 (10) 


for which the matrix At is 



( 11 ) 


Gilquin, Laurens and Rosier [1] have shown that any strictly hyperbolic 2x2 system can be 
reduced to this case by a transformation of variables, and Noelle [2] finds that any linear initial 
data is the sum of three simple waves. 

The eigenvalues of At are Aj = — l,Aa = 1, with eigenvectors 
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These can be normalised and rewritten in terms of the wave direction 9 as 


— sin 9/2 

1 

cos 9/2 

cos 9/2 

> r 2 — 

sin 9/2 


and the expansion (3) becomes 


v t 


= <*i 


sin 9/ 2 
— cos 9/ 2 


+ C*2 


cos 9/2 
sin 9/2 


This can be solved for any choice of 9:- 

ai = u t sin 9/2 - v t cos 0/2, a 2 = u t cos 9/2 + u t sin 9/2. 


( 12 ) 


(13) 


3.2.2 Elastodynamics 

The equations governing elastic waves in a uniform isotropic medium are of the form (1) with 
q = (v,p,u) T where v is the velocity, p is the trace of the strain tensor, and cj forms the 
antisymmetric (rotational) part of the strain tensor. The strains are related to the displacements 
u by 

p = div u, 


and 


u = curl u. 


The matrix A/ can be written as 


0 

0 

0 

0 

^3 

-«2 


0 

0 

0 

^2 

^3 

0 

1 1 


bi x 

bt 2 

bl 3 

0 

0 

0 

0 


0 

—(■z 

0 

0 

0 

0 


Is 

0 

-l 1 
0 
0 
0 
0 


—t 2 
tx 
0 
0 
0 
0 
0 


(14) 


Here units have been chosen so that transverse waves (S-waves) travel with unit speed, and 
longitudinal waves (P-waves) with speed y/b, where the parameter b is 

, , 2 M 

b ~ l + — 

and L , M are the Lame coefficients. 

The matrix (14) has eigenvalues Ai , 2 , 3 , 4 , 5 , 6, 7 = — — 1 — 1,0, 1,1, with non-stationary 
eigenvectors 



y/bl 


s 


r l = 

-1 

7 r 2 = 

0 

1 T3 = 


0 


-t 
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s 


r 5 = 

0 

y r 6 = 

0 

II 
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t 



t 

0 

—s 

Vbt 

1 

0 


(15) 

(16) 
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Thus the expansion (3) for this case is 


Vt 


-be 


— S 


~t 


t 


$ 


be 

Pt 

= OCX 

Vb 

+ 

0 

+ <*3 

0 

+ «5 

0 

+ <*6 

0 

+ <*7 

Vb 
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t 


0 


( 17 ) 


From this we observe that u t lies in the plane spanned by s, t , so that ( must be chosen normal 
to = curlu. For any such choice of £, 


a? 

+ 

1 

= Pt/Vb , 

a? 

- 

ax 

= l 

u t /6, 

<*6 

+ 


= t 


<*6 

- 

a 2 

= s 

Vt, 

<* 5 

+ 

<* 3 

— s 

'U t , 

<*5 

- 

<*3 

= t • 

Vt- 


(18) 


3,3 An Example with no Solution 

The following is based on a suggestion of Jonathan Goodman (Courant Institute, NYU). Con- 
sider a system comprising two copies of the acoustic equations, with unknowns (i>i,pi, V2iP2)- 
The resulting one-dimensional wave problem has a solution only for data with || 


4 Applications 

These results have potential application to the numerical analysis of wave motion, where the 
solution is often assumed to be linear within elements, and where the update strategy may 
depend on whether waves are entering or leaving some region. For a critique of methods where 
wave directions are determined relative to cell boundaries, see [3]. Dissatisfaction with such an 
approach leads to an interest in representing multidimensional waves in a way that is coordinate- 
free, but simple enough for computation. The lack of a general pattern may not be a drawback; 
exploiting behaviour that is very specific to a particular equation set may be the most effective 
approach. 

Specific applications will turn on details of numerical technique as well as the analytical 
decompositions; for that reason they are not explored here, but reference may be made to 
computations of the compresible Euler equations by splitting finite-element gradients into steady 
and unsteady components [4], and to a boundary condition for Maxwell’s equations that uses 
the above formulae to identify the direction in which waves exit the boundary [5]. 

In a non-numerical context, there is precedent for this analysis in the work of Frohn [6] who 
studied steady three-dimensional supersonic flow by the method of strained coordinates, where 
the objective is to find a coordinate perturbation such that the singularities of linear theory 
become realistic shocks and rarefactions. She selected a plane containing the streamline and the 
pressure gradient as the one in which a two-dimensional straining should be applied 

5 Conclusions 

Given a set of hyperbolic partial differential equations in D space dimensions (D > 1) and time, 
together with linearly varying initial data, one may seek directions in which the evolution of the 
problem takes place one-dimensionally. We have given examples in which the search succeeds 


5 



(and a rather artificial counterexample). Since a practical computational algorithm can only be 
based on a finite number of wave directions, this offers a possible basis for the creation of such 
algorithms. 
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